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<U , Abstract 

We compute numerically the topological charge distribution in the initial stage of a 
high energy heavy ion collision. This charge distribution is generated by Chern-Simons 
number fluctuations associated with the dynamics of strong classical fields in the initial 
state. The distribution is found to be quite narrow at RHIC and LHC energies reflecting 
a small value of the topological susceptibility. Thus the effective potential of classical 
fields is shallow in the ^-direction likely creating favorable conditions for the subsequent 
generation of P-odd bubbles. 
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1 Introduction 



In high energy heavy ion collisions, QCD matter is produced at very high energy densi- 
ties. An issue that has aroused considerable interest is the possibility [p], @ of creating 
metastable states that break the discrete global symmetries of parity (P) and charge- 
parity (CP) satisfied by the QCD Lagrangian []. There is a theorem by Vafa and Witten 
which states that there can be no stable parity violating phase of QCD [EJ. The condi- 
tions of applicability of this theorem have attracted much attention recently ||, [5], The 
Vafa- Witten theorem does not however preclude the possibility of forming metastable 
domains that violate parity (and charge-parity) 0. Furthermore, it has been suggested 
recently that the Vafa- Witten theorem may not be applicable at finite temperature ||. 

Several experimental signatures have recently been suggested that would be sensitive 
to the formation of these P and CP odd metastable states 0. Experimental searches 
are currently underway at RHIC. A preliminary report on their status was discussed at 
Quark Matter 2001 g. 

In the discussions of Refs. |], §, the formation of P and CP-odd domains is studied 
in the context of the restoration of the axial U(l) symmetry at finite temperature. A 
non-linear cr-model which incorporates the breaking of the U(l) axial symmetry Q, the Di 
Vecchia-Veneziano- Witten Lagrangian [|T0|| , (VVW) is constructed and it is shown that 
there exist metastable states (corresponding to local minima of the potential) which spon- 
taneously break P and CP The likelihood of forming long-lived metastable states 
depends on the coefficient a ~ m^, of the anomaly term in the VVW Lagrangian - they 
are more likely when a(T) — > as the temperature approaches the deconfinement tem- 
perature, T — > T c . There is evidence from lattice simulations that this is the case. The 
anomaly term is proportional to the topological susceptibility and lattice data suggest 
that the topological susceptibility for both quenched and full QCD drops rapidly in the 
vicinity of T c and is consistent with zero by T > 1.5T C (1.2T C ) for Nf = 2 (Nf = 4) flTTf . 
Further, large N c analyses suggest that a(T) — > as T — > T c [p], |], [13], [L^] . 

In Ref. @, the decay of P and CP odd metastable states into pions is described by 
the Wess-Zumino- Witten term [|14]] in the chiral Lagrangian. This term preserves the 
net chirality or handedness generated in the P and CP odd domains. The magnitude of 
the net asymmetry produced by this mechanism in heavy ion collisions is estimated to 
be ~ 10 -3 @. The signatures of P and CP odd domains were investigated also in Refs. 
T5| , |T(| |T7| ], and include an enhancement in the r\ and 7/ yields, generically associated 



with the U a\X) restoration 18 



It is of great interest to investigate how the topological susceptibility is affected by 
the dynamical conditions generated in the early stage of a heavy ion collision. Immedi- 
ately after a high energy nuclear collision, as a consequence of the phenomenon of gluon 
saturation []19 |, large color electric and magnetic fields are generated. These can in prin- 
ciple produce a large amount of topological charge through the mechanism of sphaleron 

1 This is true in the absence of a 6*-term which would break P and CP explicitly. We will assume 
that 6 — here. 

2 For a recent study of Ua(1) restoration in the linear sigma model, see Ref. 0. 
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transitions already in the initial stage of the heavy ion collision. However, as we will 
argue below, the boost invariance of the initial stage of the heavy ion collision suppresses 
sphaleron transitions. The primary mechanism for the generation of topological charge 
at the early stage is by fluctuations of the color electric and magnetic fields. For a 
translationally invariant system, the fluctuations in the topological charge determine the 
topological susceptibility. 

In this paper, we estimate the size of topological charge fluctuations in the early 
stage of a heavy ion collision. Our computation is performed within the framework of 
a classical effective field theory (EFT) |20| , approach to small x physics. In this 



approach, parton distributions in a nucleus before the collision form a color glass con- 
densate (CGC) [^, |23[. The shattering of the color glass condensate and subsequent 
multiparticle production in high energy scattering has been discussed for eA ||2U|| , pA [ 2^ 



and for nuclear collisions p4|, • For nuclear collisions, the problem can be formulated 
in the EFT |24| and the temporal evolution of gluons produced after the collisions can be 



solved for numerically PEfl . The initial energy [271] and number [2^] of produced gluons 



have been computed previously in this approach. The CGC idea has been shown to be 
consistent with many features of data from heavy ion collisions [3C, |3"Tf . 



We find in the CGC model that the fluctuations in topological charge are quite small 
at RHIC and LHC energies. The root mean squared topological charge is only one unit 
per two units of rapidity at RHIC energies and one unit per unit of rapidity at LHC 
energies. Our results suggest that the form of the effective potential in the ^-direction 
is quite shallow at early times. This result makes favorable the formation of P and CP 
odd bubbles at later times - namely, the scenario that was first discussed in Refs. ]!]. Q . 

Recently instanton/sphaleron approaches to multiparticle production have been in- 
vestigated [^2], |33], |33|, |55| . In Ref. it was shown that in the initial stage of nuclear 
collisions instanton transitions (and therefore the topological susceptibility) are severely 
suppressed by strong classical gluon fields. Shuryak has suggested that the dominant 
mechanism for particle production in nuclear collisions may be through the decay of 
sphalerons |35|| . If so, a likely consequence is that the ensuing sphaleron transitions will 
generate large P and CP violation already in the initial stage of the nuclear collision. In 
our approach, however, sphaleron transitions are suppressed by strict boost invariance. 
Sphaleron transitions, if they occur, will take place at a later stage when the strict boost 
invariance of the initial state is lost. We have estimated the topological charge generated 
by sphaleron transitions in a hot gluon plasma and find, for RHIC energies, that it is 
large compared to the topological charge generated in the initial stage by classical fields. 

The paper is organized as follows. In section 2, we briefly review the CGC picture 
of classical gluon production in heavy ion collisions. The no-go theorem for sphaleron 
transitions for a strictly boost invariant system is derived in section 3. Numerical results 
for the Chern-Simons charge are presented in section 4. Its relation to empirical observ- 
ables is also discussed there. In section 5, we discuss estimates for P and CP violation 
from sphaleron transitions at finite temperature in equilibrium hot gluodynamics. We 
also comment on the implications of these results (versus those in the boost invariant 
CGC picture) for the heavy ion experiments. 
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2 Classical gluon production in nuclear collisions 



In a high energy heavy ion collision, the dynamics of the central rapidity region is de- 
termined by the small x modes in the individual nuclei before the collision. At small x, 
the gluon density grows rapidly - however, repulsive effects soon become important and 
slow down the growth of the gluon density. This phenomenon is called saturation and 
has been discussed extensively in the literature |T9[. The occupation number of gluons 
in the saturated regime is proportional to l/a s (A s ) evaluated at a saturation scale A s . 
If A s ^> Aqcd, the occupation number is large and classical methods may be applied 



to study parton distributions in the nuclei p0| , |2"I |. It has been shown recently that 
a RG-improved generalization of this classical effective field theory (EFT) reproduces 
several key results in small- x QCD GTL 123]. The classical EFT was first applied to the 



study of collisions of large nuclei by Kovner, McLerran and Weigert J24|. The model 



as applied to nuclear collisions, may be summarized as follows. The colliding nuclei are 
idealize to travel along the light cone The high-x and the low-x modes in the nuclei are 
treated separately. Th former corresponds to valence quarks and hard sea partons and 
are considered recoilless sources of color charge For each of the large Lorentz-contracted 
nuclei (for simplicity, we will consider only collisions of identical nuclei), this results in 
a static Gaussian distribution of their color charge density p± 2 in th transverse plane 



P([p\) oc exp 



2A? 



dV t p? 2 (rt) 



The variance A s of the color charge distribution is the only dimensionful parameter of 
the model, apart from the linear size L of the nucleus. For central impact parameters, 
A s can be estimated in terms of single-nucleon structure functions f2~5fl . It is assumed, 



in addition, that the nucleus is infinitely thin in the longitudinal direction. Under this 
simplifying assumption, the resulting gauge fields are explicitly boost-invariant. 

The small x fields are then described by the classical Yang-Mills equations 

D^ v = J u (1) 

with the random sources on the two light cones: J u = J2i2^,±^{ x ^)Pi,2( r t)- The two 
signs correspond to two possible directions of motion along the beam axis z. As shown 
by Kovner, McLerran and Weigert (KMW) f24|, low-x fields in the central region of the 
collision obey sourceless Yang-Mills equations (this region is in the forward light cone of 
both nuclei) with the initial conditions in the A T = gauge given by 

A* = A\ + A\- A ± = ±^x ± [A\, A\\. (2) 

Here the pure gauge fields A\ 2 are solutions of flU) for each of the two nuclei in the 
absence of the other nucleus. 

In order to obtain the resulting gluon field configuration at late proper times, one 
needs to solve (ffl) with the initial condition (□). Since the latter depends on the random 
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color source, averages over realizations of the source must be performed. KMW showed 
that in perturbation theory the gluon number distribution by transverse momentum (per 
unit rapidity) suffers from an infrared divergence and argued that the distribution must 
have the form 

1 M.,\ 4 . (k\ 



for k± ^> A s . The log term clearly indicates that the perturbative description breaks 
down for k± ~ A s . 

A reliable way to go beyond perturbation theory is to re-formulate the EFT on 
a lattice by discretizing the transverse plane. The resulting lattice theory can then be 
solved numerically. We shall not dwell here on the details of the lattice formulation, which 
is described in detail in Ref. |26|, |27||. Keeping in mind that A s and the linear size L of the 
nucleus are the only physically interesting dimensional parameters of the model [ 22] , we 
can write any dimensional quantity q as Af/ g (A s L), where d is the dimension of q. All the 
non-trivial physical information is contained in the dimensionless function f q (A s L). We 
can estimate the values of the product A S L which correspond to key collider experiments. 
Assuming Au-Au collisions, we take L = 11.6 fm (for a square nucleus!) and estimate 
the saturation scale A s to be 2 GeV for RHIC and 4 GeV for LHC ||25|| . Also, we have 
approximately g = 2 for energies of interest. The rough estimate is then A S L ph 120-150 
for RHIC and A S L ph 240-300 for LHC. Since the gluon distribution in nuclei is not known 
to great precision, there is a considerable systematic uncertainty in these estimates. 

This model has been applied recently to study classical gluon production. The energy 
and number distributions have been computed numerically, and the dependence of these 
quantities on A s has been determined |27|, ^8] . We will now apply the model to compute 
the Chern Simons number generated in the early stages of a heavy ion collision. 



3 Boost invariance and Chern Simons number 



In general, the Chern-Simons (CS) number changes under a gauge transformation by 
an integer equal to the winding number of that transformation. Since the potential of 
the theory is invariant under all gauge transformations, there exists a direction in the 
functional space along which the potential is periodic. The system can move along that 
direction arbitrarily far from its initial configuration. Such motion generates arbitrarily 
large changes of the CS number or, equivalently, arbitrarily large values of topological 
charge. Dynamical, real time, processes corresponding to nontrivial gauge transforma- 
tions are called sphaleron transitions. In a large space-time volume these transitions 
proceed independently in causally unrelated regions. If, in addition, the space-time is 
translation- invariant, the topological charge for a large space-time volume is a random 
walk: 

((N cs (t)-N cs (t')) 2 )=TV(t-t>), (4) 
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where T is known as the sphaleron transition rate, or as the topological susceptibility. 
Here, and in the following, (...) means an ensemble average: the thermal ensemble in 
thermal equilibrium, and the ensemble of (central) nuclear collision events in the case of 
interest here. This random walk behavior was indeed observed in numerical simulations 
of SU(2) and SU(3) Yang-Mills theories at high temperature. 

The dynamics of the CS number in the central region of a nuclear collision differs from 
the generic case in two ways. Firstly, the time translation invariance is broken by the 
instance of the collision itself. Secondly, and more importantly, there are no nontrivial 
boost-invariant gauge transformations. Consequently, for boost-invariant configurations, 
such as those in the central region, the potential is not periodic, and there can be no 
sphaleron transitions. Instead, one expects the CS number to fluctuate in the vicinity of 
zero. 

In order to make this second point explicit, we shall now examine the form of the CS 
number functional for boost-invariant configurations, and show it to be invariant under 
all boost-invariant gauge transformations. It will then follow that, in this special case, 
the CS number cannot be changed by a dynamical process equivalent to a gauge trans- 
formation, i.e., that sphaleron transitions are prohibited for a boost invariant system. 

It is convenient for us to use the coordinates (f t ,T,rj). Here r = a/2x + x~ is the 
proper time, t] = ln(x + /x~)/2 is the space-time rapidity, and f t = (x,y) is the position 
in the transverse plane. For the components of the gauge potential, the boost invariance 
assumption (see the discussion immediately prior to Eq. |l|) means that 

A t (r, x t , rj) = A t (r, x t ) ; A v (t, x t , rj) = $(r, x t ) , (5) 

and the system effectively becomes (2+l)-dimensional, with A v becoming an adjoint 
scalar field. It is also convenient to remove the fourth component of A by imposing the 
gauge condition A T = 0. 

Now consider the quantity 



l — J d 2 x t $ a B a , (6) 



16-7T 

where 

B a = F a xy = d x A; - d y A a x + ge ahc A\A% (7) 

is the component of the magnetic field perpendicular to the transverse plane. This 
quantity is obviously invariant under all rapidity independent gauge transformations. In 
order to verify that v is the CS number per unit rapidity, we will compute its proper time 
derivative (denoted by a dot in the following) and show it to be equal to the topological 
charge density integrated over the transverse plane. We recall that the topological charge 
density is given by 

-^E-B, (8) 
16tt 2 ' V ; 
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E a and B a being the color electric and magnetic fields, respectively; the anomaly equation 
relates this quantity to the divergence of the chiral current. 

To this end, it suffices to write out v explicitly, using (|7|), and integrate the terms con- 
taining d x Ay and d y A® by parts in the transverse plane, assuming that the corresponding 
boundary terms vanish. The result is 

where is the covariant derivative: D^Q = di<§> a -\-e A\<& c . To complete the proof, we 
note that (a) for boost invariant configurations, the integrand of (0) is equal to rE a ■ B a ; 
and (b) the space-time volume element in the r,7],r t coordinates is d 4 x = rdr/drd 2 ^. 
Hence we have shown that the proper time variation of v in a boost invariant process gives 
the topological charge per unit rapidity of the process. If any such process is equivalent 
to a gauge transformation, it leaves v unchanged, the corresponding topological charge 
vanishes, and no sphaleron transitions are possible. 



In the EFT description of nuclear collisions the scalar field $ vanishes at r = p6 
and, therefore, v{j — 0) = 0. In other words, the topological charge of the process at 
any given r is simply tji/(t), where rj is the rapidity extent of the process. Since the 
EFT has a parity-even action, both signs of v(t) are equally probable, and the ensemble 
average of u(t) vanishes. Our principal object of interest is then (y 2 (r)). In the absence 
of sphaleron transitions, we cannot expect the proper time evolution of v to resemble a 
random walk. Indeed, as we shall see in the next section, (^ 2 (t)) approaches a constant 
at large r. 

Strictly speaking, because of the strict boost invariance built into our approach, 
(z/ 2 (r)) per unit space-time volume cannot be interpreted as the topological susceptibility. 
Nevertheless, the computed (u 2 (t)) is a measure of the contribution of classical fields to 
the topological susceptibility. As will be discussed further in Section 5, the likelihood of 
forming long-lived P and CP odd domains depends sensitively on this quantity. 



4 Numerical results 



We use a lattice formulation of the EFT developed in |26|, [27], |28fl . As explained in these 
references, the key dimensionless parameter for those observables which exist in the 
continuum limit is A S L, where L is the radius of the nucleus (assuming a square nucleus, 
L 2 = ttR 2 ). The value of A s , expressed in units of the lattice spacing a, indicates 
how close the system is to the continuum limit. In the current study, for an SU(2) 
gauge theory Q we set A s a = 0.29. The CS number v is a dimensionless quantity, whose 

3 We study an SU(2) gauge theory since it is simpler to simulate numerically. The physically more 
interesting case of an SU(3) gauge theory was studied recently and gluon number distributions com- 
puted No significant qualitative differences from the SU(2) case were found for the energies of 
interest. 
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A s r 

Figure 1: Average squared topological charge vs the proper time (in units of A^ 1 ) for 
A S L = 74.2 (diamonds), 148.4 (pluses), and 297 (squares). The latter two values of A S L 
correspond to the RHIC (A s = 2 GeV) and LHC (A s = 4 GeV) regimes respectively. 
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Figure 2: Probability distribution of v in the proper time interval 30 < A s r < 65 for 
A S L = 148.4. These corresponds to 3 fm < r < 6 fm for an estimated value A s = 2 GeV 
for RHIC energies. 
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ultraviolet properties are similar to those of the total gluon number N, studied in [28[] . In 
that study, we determined that discretization artifacts are very small for N at A s a = 0.29. 

We examined the proper time evolution of v for three values of A S L, namely, 74.2, 
148.4, and 297. The last two of these correspond roughly to central Au-Au collisions 
in the RHIC and LHC regimes, respectively. In each of the three cases, we generated 
about 80 independent initial configurations and solved the classical equations of motion 
for the duration of proper time of at least 60/A s (this corresponds to 6 fm for RHIC). 
The average proper time histories of v 2 in each case are presented in Figure 1. Evidently, 
v varies over a much longer time scale, compared to other global quantities, such as the 
total energy and the total particle number |27|, p8"| . Nevertheless, (z/ 2 (r)) does approach 
a constant value at late r. That constant value is a rapidly growing function of A S L. 
Since the dynamics considered here is exactly boost-invariant, ((^ 2 (t))) 1//2 is a measure 
of the anomalous chirality violation per unit rapidity in the initial stage of a central 
collision. 

To estimate the magnitude of the effect, we took the proper time average of (z/ 2 (r)) 
between A s t = 30 and A s t = 65. The results are summarized in Table 1. If we 
take the experimental range of two units of rapidity, it follows that for RHIC we have 
approximately one unit of the topological charge squared for two units of rapidity while at 
LHC one may have two units of topological charge squared in the same rapidity interval. 
Note that (u 2 ) is approximately proportional to (A S R) 2 , i.e., to the transverse area in 
units of the nonperturbative scale A s . This means that dynamical processes contributing 
to the topological charge become uncorrelated at spatial separations exceeding A s , and 
the resulting fluctuations of v add up randomly. 



A S R 


41.9 


83.7 


167.5 


(u 2 ) 


0.046 ± 0.004 


0.16 ±0.01 


0.69 ±0.05 



Table 1: The average value of (v 2 ), tabulated as a function of the dimensionless parameter 
A. S R. 



As we have explained in the preceding sections, there is no CS number diffusion in 
a boost invariant setting, nor is there a ground state degeneracy. It is for this reason 
that (z /2 (t)) approaches a constant at large r, instead of growing linearly, as would be 
the case if sphaleron transitions were permitted. Also, as Figure 2 illustrates, there is 
no special significance to integer values of v other than zero, which in a general case 
correspond to degenerate minima of the potential. 
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5 Comparison with other estimates of the topologi- 
cal susceptibility and experimental ramifications 



Thus far, we have assumed that strict boost invariance holds in heavy ion collisions 
at very high energies. If strict boost invariance does not hold, sphaleron transitions 
are allowed. Whether sphaleron transitions turn on slowly or rapidly when strict boost 
invariance is relaxed is not clear. A way to answer this question is to introduce a small 
perturbation a(x t ,r),T), namely, 

A(x t , r], r) = A(x t , t) + a(x t , 77, r) , (10) 

in the Yang-Mills equations and to study how it evolves. While the Yang-Mills equations 
equations can be linearized with respect to a(x t , 77, r), such a study has to be performed 
numerically, since A(x t ,r) is not known analytically [jSJ. 

Let us assume that the system does thermalize and strict boost invariance of field 
configurations is lost []. We then need an estimate of the sphaleron transition rate in 
a hot gluon plasma. To simplify the setup, we shall first ignore dynamical fermions 
and assume that the temperature T is much larger than any other dimensional scale 
in the problem. It then follows from dimensional considerations |38] that the rate is 



proportional to T 4 However, the dependence of the rate on the coupling constant g is 
less obvious. In the weak-coupling regime the sphaleron transitions are dominated by 
soft modes, with momenta of the order of g 2 T. The time scale for the transition depends 
on how important is the interaction between these degrees of freedom and hard thermal 



modes. A series of perturbative studies [44, Hq], performed at extremely weak coupling 



(ln(l/g) ^> 1), suggest that this interaction sets the time scale for sphaleron transitions 
to ((7 4 T) _1 up to logarithmic corrections. Then the coupling constant dependence is 
T oc ajjT 4 , and a recent numerical estimate gives [f37j 



sph 



108 a, T fir 



» 



On the other hand, there are indications that, away from the asymptotically weak cou- 
pling regime, the relevant time scale for the soft modes is (<7 2 T)" 1 [|46 |. Should this be 



the case, the numerical estimate [37] should be re-interpreted to give 



r sp ,^8« 4 T 4 . (12) 

In the following, we will use both estimates. 

For a hot plasma (without dynamical fermions) the potential is periodic and one 
generates topological charge via a random walk. One has 

(AQjj) = [ T d 4 xT sph , (13) 
J 



4 Boost invariance may still be preserved when one averages over events. 



10 



where d 4 x = d 2 Xtrdrdr]. In an expanding boost invariant plasma, the temperature 
decreases with time. For an isentropic expansion one has T oc r -1 / 3 . More precisely, 



90 dS 1 1 



4-7T 2 dr\ L 2 gd r 



1/3 



(14) 



where g& is the degeneracy factor and dS/dr] is the entropy per unit rapidity. From 
Eq. O and Eq. IT3L we obtain 



drj 



8 



dS 1 90 



d-q g d 47r 2 



4/3 



(15) 



where A and n are correspond to the different coefficient and power respectively in 
the two sphaleron transition rates. At RHIC energies, dS/drj ~ 3.6 x 1000, a s ~ 0.3, 
Tf ~ L ~ 10 fm, gd = 16. Then from the sphaleron transition rate in Eq. [TTj, and from 
Eq. ITR we find for an expanding plasma, 



rhic » 1600 ; ~ " Z°' RMS \rhic ~ 40 . (16) 



d(Ag§)W ~ 1600 . d < A <?5>gjl 

j \rhic ~ iouu , — 

ct?7 a?7 

Here RMS denotes root mean square. Similarly, for the sphaleron transition rate in 



Eq. [12|, and Eq. [15, we have 



; 1 RHIC ~ 4 UU , |_R/f/c~^U. l^l/J 

0(77 a?7 

Though these estimates of the topological charge generated by sphaleron transitions differ 
from each other they both are significantly larger than the topological charge squared 
obtained from boost invariant classical fields at RHIC (and LHC). In this paper, we 
have considered the case of pure gluodynamics; the effects due to dynamical fermions 
had been addressed in Refs. g0|, §]], @. 

We now come to the phenomenological implications of our results. As discussed 
previously, the likelihood of long-lived P and CP violating metastable states in heavy 
ion collisions [l], ^[ depends sensitively on the coefficient a of the anomaly term in the 
WW Lagrangian. This coefficient is proportional to the topological susceptibility. The 
likelihood of metastable states is greater for smaller values of a. Lattice and large N c 
calculations suggest that the topological susceptibility is indeed small for T > T c The 
previous estimates did not include the possibility of contributions from fluctuations in 
classical fields or from sphaleron transitions. A difficult problem is how to convert 
the estimates for generating non-zero topological charge into one for P and CP odd 
asymmetries for measured pions |. We will not address this problem directly here. 
Instead, since the quantity (v 2 ) that we compute here is closely related to the topological 
susceptibility, we can relate our estimate to those of Refs. |], ||. Our estimates here 

5 We speculate that one may construct, akin to L • S coupling, invariants resulting from the coupling 
of the vorticity of CS number to that of the vorticity of collective flow in an expanding plasma. 
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suggest that the contribution to a from classical fields is small and would not therefore 
affect the estimate of Ref. ||. On the other hand, the contribution from sphaleron 
transitions might be potentially large. Thus if the RHIC experiments see a large P and 



CP violating effect, it may also be due to sphalerons |35 



Our results depend on the strong boost invariance assumption, namely, each event 



is boost invariant. At central rapidities, the measured spectra look boost invariant |49| 
on average. However, what is relevant to our assumption is the requirement that the 
spectra are boost invariant on an event by event basis. One way to ensure this is to 
compute the standard deviation from strict boost invariance in one unit of rapidity 
about mid-rapidity and check if it is small. Such studies are underway but no conclusive 



results have been published yet ||50|| . As one goes to larger rapidities, the boost invariance 
assumption breaks down. The no-go theorem forbidding sphaleron transitions is violated 
and one may therefore see larger effects at non-central rapidities (in the P and CP odd 
observables currently being investigated) due to sphaleron transitions. It will be very 
interesting to see if this prediction can be tested in the on-going experiments at RHIC. 
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